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Bulk and surface plane electromagnetic waves in 

anisotropic media 

Abstract 

A new analytical approach to description of electromagnetic waves in nonmagnetic 
anisotropic media is presented. Amplitudes of their reflection and refraction at interfaces 
and also reflection and transmission of plane parallel plates are derived. Beam splitting 
at reflection, and creation of surface waves at the interfaces are studied. A simple labo- 
ratory demonstration of the beam splitting is proposed. D'yakonov surface waves, their 
description and observation are discussed. 

— Reading the manuscript, this reviewer decided that 
the authors must be new to electromagnetics. Oth- 
erwise, they would not have wasted their time rein- 
venting the wheel. Everything that could be known 
about propagation in an uniaxial material has been 
known since at least 184.1- So there is no sense wast- 
ing time there. 

Referee of Opt.Technol. 



1 Introduction 

Description of electromagnetic waves in homogeneous anisotropic media did not change since 
Fresnel times more than 160 years ago. Here the first time since then we present a new approach, 
which makes physics here very simple and transparent. We will not waist time describing the 
standard approach, which can be found with only slight variation in all the textbooks on 
electrodynamics or optics [1-13] containing chapters on anisotropic media. Instead we directly 
start with our approach. 

An anisotropic medium is characterized by some direction, called axis which we will denote 
by a unit vector a. In a plane electromagnetic wave 

£ exp(ifcr — iut) (1) 

propagating in such a medium at an arbitrary direction with respect to a the wave vector k{u!) 
and polarization vector £ depend on angle 9 between a and direction of propagation k — k/k. 

The main feature of our approach (similar to the one used for elastic waves [14]) is a special 
representation of the dielectric permittivity tensor e. It was proven by Fedorov [12] that in a 
uniaxial anisotropic medium the tensor e can be represented as 

£ij = eiSij + e'aiUj, (2) 

where Ci is isotropic part, and anisotropy is characterized by the unit vector a with components 
Gi and by anisotropy parameter e'. 

In the next section we find £ and k{u)) in anisotropic media with one and two mutually 
orthogonal axes. In the third section we discuss refiection of these waves from interfaces between 
anisotropic and isotropic media, study beam splitting at refiection, conditions for creation of 



surface waves, and propose a device to demonstrate the beam splitting at reflections in a 
laboratory. In section 4 we calculate reflection and transmission for plane parallel anisotropic 
plates. In 5-th section we calculate speed of the D'yakonov surface waves [18], which are an 
analog of Rayleigh elastic waves on a free surface of elastic media. We correct some defect of 
the derivation of these waves in [18] and propose an experiment to generate and to observe the 
D'yakonov waves. In 6-th section we summarize our results. In sections 3-5 we limit ourselves 
only to uniaxial media and only explain the ideas. An approach to biaxial media and most 
heavy mathematics are shifted to appendices. 



2 Plane waves in anisotropic media 

— / don't get the goal of this paper 
Referee of Am.J.Phys. 

The wave equation is derived from Maxwell equations, which in the absence of currents and 
charges are 

-[VxE{r,t)] = —B{r,t), [V x H{r,t)] = —D{r,t), VB = 0, VD = 0, (3) 
where 

B^fxH, D = eE, (4) 

and II, £ are magnetic and dielectric permittivities. In the following we take /i — 1, and then 
the equations (3) are simplifled to 

-[VxE{r,t)]^—H{r,t), [V x H{r,t)] ^ —£E{r,t), VH^O, VsE^O. (5) 

Differentiation of the second equation of (5) over time and after that substitution of the first 
equation leads to the wave equation 

-[V X [V X Eir, t)]] = J^,eE{r, t). (6) 
2.1 Electric vector of the wave 

Substitution of the plane wave (1) into this equation reduces (6) to 

k'^S - k{k ■ S) = kleS, (7) 

where ko = uj/c. Eq. (7) is valid only in homogeneous media. If we have an interface between 
two homogeneous media we have two different equations of the type (2) in them, and matching 
of waves in two media is performed via boundary conditions, which follow not from the wave 
equation itself, like in quantum mechanics, but from Maxwell equations. 

In uniaxial anisotropic media we use tensor e in the form (2). Therefore for a plane wave 
(1) we have 

e£ = ei£ + e'a(a-£), (8) 
and the last equation in (5) leads to 

ei(fc-£) + e'(fc-a)(a-£) = 0. (9) 
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Substitution of (8) into (7) gives 

k^£ - k{k ■ S) - kleS = {k^ - klei)S - k{k ■ £) - kle'a{a ■£)=0. (10) 

To find £ we need to solve (10) with account of (9). 

The 3-dimensional vector £ can be represented by coordinates in some basis. If k is not 
parallel to a, we can use as a basis three independent vectors a, k — k/k and 

Ci = [a X k\. (11) 

In this basis £ looks 

£ ^ aa + pK, + (12) 

with coordinates a, (3 and 7, which are not independent, because of Eq. (9). 
Substitution of (12) into (9) gives 

ei [k^ + {k-a)a]+ e\k • a) [a + /3(k • a)\ = 0, (13) 

from which it follows that 

l + ri[K-ay 

where 77 = e'/ei. Substitution of (14) into (12) gives 

/ (K-a)(l + r7)\ , , 

£^a\a-K, \ ^ +7ei = Q;e2 + 7ei, 15 

which shows that £ hes in a plane of two independent vectors ei = [a x n] and the orthogonal 
to it 

1+77 

62 = a - k(k • a)— — = a - k(k • a)e2{9)/ei, (16) 

1 + 77(k ■ ay 

where cos 6 = k ■ a and we introduced anisotropic dielectric permittivity 

62(g) -6. + (17) 
1 + 77 COS"' 9 

To find coordinates a and ^ we substitute (15) into (10) and multiply it by e\. As a result 
we obtain 

{e - kler)^ei = 0. (18) 
It shows that if 7 7^ 0, then (18) can be satisfied only when 

k'^ = kle^. (19) 
Multiplying (10) by a and taking into account that 

l-{K-af l-{K-af 
" • = l + r^in-ar " " = " " ^ + 7,(/.-a)^ = "^^^ " " ^'^^ ^''^ 
we obtain 

[k^ - kle2{9)) a{a ■ 62) = 0. (21) 

Therefore, if a 7^ 0, and k, Eq. (21) can be satisfied only, when 

k''^kle2{9), (22) 

where e2{9) is given in (17). Since the length of k is different for two polarization vectors, 
therefore a single plain wave can exist only with a single polarization along either 62 or ei. 

We will call "transverse" the mode with polarization £1 = Ci, and "mixed" the mode 
with polarization along £2 = 62- The mixed mode according to (20) contains a longitudinal 
component along the wave vector k. We think that such a nomenclature is better than common 
names: "ordinary" for £1 — ei, and "extraordinary" for £2 — e2, because our names point to 
physical features of these waves. 
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2.2 Magnetic fields 

Every electromagnetic wave besides electric contains also magnetic field. From the equation 
V • If = 0, which is equivalent to k ■ H — 0, it follows that the field H is orthogonal to fc. It 
is also orthogonal to £, which follows from the first equation of (5). After substitution of Eq. 
(1) into (5) and the field H in the plane wave form 

H{r,t) ^nex.p{ik-r -iut), (23) 

with polarization vector 7i we obtain 

n^^[Kx£]. (24) 

Kq 

For transverse and mixed modes in uniaxial media, respectively, we therefore obtain 

Hi = —[k X ei] = —[k X [a X k]], H2 = t^[k x 62] = x a], (25) 

Kq Kq Kq Kq 

and the total plain wave field looks 

*(r", t) = ipj exp{ikj ■ r — iut), (26) 

where = + and j denotes mode 1 or 2. In isotropic media we also can choose, say 
£ = [a X k] and 'H = [k x [a x k]]. However, there a can have arbitrary direction, therefore 
the couple of orthogonal vectors £ and H can be rotated any angle around the wave vector k. 

3 Reflection from an interface between uniaxial and 
isotropic medium 

— / spent a few separate sittings reading this paper 
hut cannot convince myself to sit with pencil and pa- 
per and follow the mathematics presented here, since 
I know I can open Jackson,... or Griffith and get 
really all I would need to know about the interac- 
tion of E&M waves with the surface of a dielectric 
medium-in a more streamlined, concise and easy to 
follow manner. 
Referee of Am.J.Phys. 

A note Anisotropic media are not considered in [15,16]. 

Imagine that our space is split into two half spaces. The part at 2; < is a uniaxial 
anisotropic medium, and the part at 2; > is vacuum with ei = 1, 77 = 0. We have two 
different wave equations in these parts, and waves go from the reign of one equation into the 
reign of another one through the interface where they must obey boundary conditions imposed 
by Maxwell equations. 

Let's look for reflection of the two possible modes incident onto the interface from within 
the anisotropic medium. 

3.1 Nonspecularity and mode transformation at the interface 

First we note that reflection of the mixed mode is not specular. Indeed, since direction of 
k after reflection changes, therefore the angle 9 between a and k does also change, and k, 
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according to (22), changes too. However the component k\\ parallel to the interface does not 
change, so the change of k means the change of the normal component k±, and this leads to 
nonspecularity of the reflection. 

Let's calculate the change of k± for the incident mixed mode with wave vector k2r, where 
the index r means that the mode 2 propagates to the right toward the interface. For a given 
angle 9 between k2r and a we can write 



k 



2r± 



\l + r]cos^e ^11' ^^^^ 



however the value of k2r± enters implicitly into cos^, so to find explicit dependence of k2r± on 
a it is necessary to solve the equation 

fcj + + rjikiiH ■ a) + x{n ■ a)f = A;^ei(l + rj), (28) 

where x denotes k2r±, n is a unit vector of normal, directed toward isotropic medium, and / is 
a unit vector along fey, which together with n constitutes the plane of incidence. Solution of 
this equation is 

-rjk\i{n ■a){l-a) + ^e^kl{l + r/)(l + 7]{n ■ af) - A;J(1 + r]{l ■ af + r]{n ■ a^) 
^ " l + r]{n-af " ^^^^ 

The sig n chosen b efore square root provides the correct asymptotics at 77 = equal to isotropic 
value sjeikl — /cy . 

In general vector a is representable as a = an + /3Z + 7^, where t = [nl] is a unit vector 
perpendicular to the plane of incidence. The normal component k2r± depends only on part of 
this vector a' — an + which lies in the incidence plane. If we denote a = \a'\ cos{9a), (3 = 
I a' I sin(^a), where \a'\ is projection of a on the incidence plane, and introduce new parameter 
rj' — r}\a'\'^ < rj, then formula (29) is simphfied to 



sin(2^„) + 2^eiki{l + r])[l + r]' cos2(^„)] - A;J(1 + f]') 

For the reflected mixed mode (mode 2, propagating to the left from the interface) an equation 
similar to (28) looks 

kl + x'^ + r]{kii{l ■ a) - x{n ■ a)f = A;gei(l + 77), (31) 
where x — k2i±, and its solution is 

Tj'k^i sin(2^„) + 2^eiA;2(l + v)[^ + V' cos^{Oa)] - ^J(l + V') 



^2l± 



2[l + ?7'cos2(^J] 



(32) 



We see that the difference of the normal components of reflected and incident waves of mixed 
modes k2i± — k2r± is 

rj'k\\ sm{29a) 

k2i± - k2r± = ^7^. (33) 

1 + 77' COS'' [9a) 

In the following we will present such differences in dimensionless variables 

_ k2l± - k2r± ^ V'QM'^Oq) + V(l + + V' COS^(^a)] " Q'{^ + V') 

koV^i 2[l + 77'cos2(^„)] ' ^ > 



-0.4 
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Figure 1: Variation of change A of normal components for reflected and incident waves in 
dependence of angle 9 = 9a of anisotropy vector with respect to normal n. Vector a is supposed 
to lie completely in the incidence plane. The curves Aij represent dimensionless ratio A.ij{9) 
given by (34), (35) and (36). The curves were calculated for rj — rj' — 0.4 and q — k\\/kQy/e{ — 
0.7. 



where (f = kyklei. The reflection angle depends on orientation of anisotropy vector a and it 
can be both larger than the specular one, when 9a > 0, or smaller, when 9a < 0. 

In the case of transverse incident mode the length k — \k\ of the wave vector, according to 
(19), does not depend on orientation of a, therefore this wave is reflected specularly. 

Every incident mode after reflection creates another one, because without another mode it 
is impossible to satisfy the boundary conditions. Let's look what will be the normal component 
of the wave vector of other mode. If the incident is the wave of mode 2, reflected transverse 
mode (mode 1 propagating to the left, away from the interface) will have kii±_ — yjeik^ — . 
Therefore according to (30) the difference A12 = {kii± — k2rA_) /ko^/e[ is 



^12 



Vgsin(2^„) + 2^(1 + r])[l + Vcos2(^„)] -q^{l + rj') 
2[l + 7;'cos2(^„)] 



(35) 



In the opposite case, when the incident mode is transverse one, the reflected mixed mode will 
have k2i± shown in (32). Therefore the difference A21 = {k2i± — kir±) / ko-y/ii is 



^21 



rj'q sm{29a) + 2^(1 + rj) [1 + y' cos^(g„)] -qH^+V) r— 
2[l + r)'cos^(9a)] ^ 



(36) 



The changes of normal components with variation of 9a according to (34), (35) and (36) for 
some values of dimensionless parameters rj and q and vector a lying completely in the incidence 
plane, are shown in Fig. 1. From this flgure it is seen that the strongest deviation of reflected 
wave from specular direction is observed for reflection of mixed to mixed mode. 

Since reflection of mode 2 is in general nonspecular, it can happen that the wave vectors of 
reflected and transmitted waves will be arranged as shown in flg. 2, and it follows that there are 
two critical angles for 93. The flrst critical angle (fd {q^ = 1/ei) is the angle of total reflection. 
The transmitted wave at it becomes evanescent. The totally reflected fleld contains two modes. 
At the second critical angle </7c2, when q is in the range 



^ ^ 2 ^ (l + 77)(l + 77W(g„)) 

1 + 77' 



(37) 




Figure 2: Arrangement of wave vectors of all the modes created by the incident wave of mode 
2, i.e. of polarization when the anisotropy vector a has the direction shown here. The 
grazing angle of the reflected mode 2, £ 2, is less than specular one (specular direction is shown 
by the broken arrow), and the grazing angle ipi of the reflected mode 1, £ 1, is even lower. The 
grazing angle (/Pq of the transmitted wave ^0 is even lower than Lpi. Wc can imagine than at 
some critical value = Lpd the angle ipo becomes zero. It means that at (/? < Lpd transmitted 
wave becomes evanescent and all the incident energy is totally reflected in the form of two 
modes. More over, there exists a second critical angle (pc2, when </?i = 0. Below this angle at 
(p < (pc2 the mode £ 1 also becomes evanescent. In this case all the incident energy is totally 
reflected nonspecularly in the form of the mode 2. At the same time the two evanescent waves 
"^0 and ^2 combine into a surface wave, propagating along the interface. The arrows over £ 
show direction of waves propagation with respect to the interface. In the figure there is also 
shown the basis which is used along the paper. It consists of unit normal vector n along normal 
(2;-axis), unit tangential vector I (x-axis) which together with n defines the incidence plane, 
and the vector t (y-axis) looking toward the reader, which is normal to the incidence plane. 




Figure 3: Dependence of dimensionless normal components of incident and refiected waves 
on g = kcosif/ko^fei. The solid curve corresponds to the incident wave moving to the right 
kr{q) = k2r±/kQy/ei. The dotted curve corresponds to the reflected wave of mode 2 moving 
to the left kl{q) = A;2/±/^o\/^- -^^^ broken curve corresponds to the reflected wave of 
mode 1 moving to the left kl{q) = kii±/koy/e{. It is seen that at g > 1 the mode 1 ceases to 
propagate. Its normal component kl{q) — ^\ — q^ — —i^q^ — 1 becomes imaginary, therefore 
the reflected mode 1 becomes an evanescent wave. Together with transmitted wave, which 
becomes evanescent at g^ = 1/ei, the mode 1 constitute the surface electromagnetic wave. 
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the reflected mode 1 also becomes evanescent. Together with evanescent transmitted wave 
the mode 1 constitutes a surface wave, propagating along the interface. In that case we have 
nonspecular total reflection of the single mode S2 and the surface wave tied to it. 

In figure 3 it is shown how do the normal components of wave vectors change with increase 
of g, which is equivalent to decrease of </?. For e\ — 1.6 the first critical angle corresponds to 
q 0.8. The second critical angle corresponds to g = 1. 



3.2 Demonstration of the beam splitting with the help of a birefrin- 
gent cone 

— In order to ensure that our journal continues to 

publish articles that cater for our broad readership, 
every paper submitted must meet our stringent edi- 
torial criteria. We believe that your article does not 
meet these criteria, so it has been withdrawn from 
consideration. 
Editor of J.Phys.A Optics 




Figure 4: Demonstration of the beam sphtting of light in a birefringent cone. Bright spots on 
a vertical screen change their position and brightness when the cone is rotated around vertical 
axis. 

Beam splitting at interfaces of an anisotropic medium can be spectacularly demonstrated 
with the help of birefringent cone as shown in fig. 4. In the geometrical optics approximation 
a narrow incident beam of light after refraction on the side surface of the cone is split into 
two rays of two different modes 1 and 2. Both modes are further split into two components at 
reflection from basement of the cone. Four resulting beams after refraction at the side surface 
go out of the cone and produce on a vertical screen four bright spots. Their positions and 
brightness depend on direction of the anisotropy axis inside the cone and vary with the cone 
rotation. 

The direct numerical calculations for parameters ei = 1.6, rj = 0.8, a in the figure plane, 
and sincK = 0.5, sin^ = 0.3 sin 7 = 0.5 show that outgoing beams from below to top have 
directions characterized by tg^i = 0.2, tg52 = 0.4, tgSs — 0.6 and tgS^ — 0.7 respectively. 



4 Calculation of refraction at interfaces and scattering 
on plates 

— Too much of a mathematical exercise and too lit- 
tle physics for Physica Scripta. 
Referee of Phys. Scripta 

To calculate reflection and transmission of a plane parallel anisotropic plate, placed in 
isotropic (for instance, vacuum) medium, it is necessary to know reflections and refractions at 



8 



interfaces from inside and outside anisotropic medium, which is obtained by imposing boundary 
conditions stemmed from Maxwell equations. Knowledge of everything at interfaces permits to 
write directly reflection and transmission of the plate by the method which is explained in [17] 
and will be shortly described below. 

4.1 Reflection and refraction from inside of the anisotropic medium 

The wave function in the full space is 

where © is a step function equal to unity, when inequality in its argument is satisfied, and to zero 
in opposite case, half space 2; < is occupied by anisotropic medium, and the half space z > 
is vacuum, = £ + H, arrows show direction of waves propagation, ip • denotes the incident 

wave of mode j (j = 1,2), ip j, (/ = 1,2) denotes reflected wave of mode j', k j = (fcy, /c^v^), 
k j' — {k\\, —kjfi±), ko — (fcy, yJkQ — /cy), ipe,m, 'T^e,mj Sire the refracted fields and refraction 
ampfitudes of TE- and TM-modes respectively for the incident j-mode. To find refiection 
~^ and refraction 1^ amplitudes (the arrow over them shows the direction of propagation of 
the incident wave toward the interface), we need to impose on (38) the following boundary 
conditions. 

4.2 General equations from boundary conditions 

Every incident wave field can be decomposed at the interface into TE- and TM-modes. In TE- 
mode electric field is perpendicular to the incidence plane, £ oat, therefore contribution of j-th 
mode into TE-mode is [Sj ■ t). In TM-mode magnetic field is perpendicular to the incidence 
plane, "H oc t, therefore contribution of j-th mode into TM-mode is {Hj ■ t). For refracted field 
in TE-mode we accept = t, — [1^0 x and for refracted field in TM-mode we accept 

"Hm — t, '^m — ~[>^0 X t\- 

4.2.1 TE-boundary conditions 

In TE-mode for incident j-mode we have the following three equations from boundary condi- 
tions: 

1. continuity of electric field 

(t ■ t,) + {t ■ ti)tlj + {t ■ t2)t2j = {t ■ te)l^ej, (39) 

2. continuity of magnetic field parallel to the interface 

(Z • %) + (/ • ^l)^i,- + (/ • %)t2j = {I ■ K X t])l^ej = -/^0±^ej, (40) 

3. and continuity of the normal component of magnetic induction, which for /x = 1 looks 

(n • itj) + (n • ^l)7"lj + (»^ • ^2)~^2j = (»^ • [f^O X t])l^ej = «0||Tej- (41) 

The last eq. (41) is, in fact, not needed, because it coincides with (39). 



9 



4.2.2 TM-boundary conditions 

In TM-modc wc have the equations 

1. continuity of magnetic field 

{t ■ %) + {t ■ + {t ■ = -t^,, (42) 

2. continuity of electric field parallel to the interface 

{1 ■ t,) + {I ■ ti)ty + {I ■ t2)t2j = -{I ■ K X t])^^j = Ko±l^mj, (43) 

3. and continuity of the normal component of field D 

(n • e'ij) + (n • £^i)~^ij + (n • £^2)~^2j = • [«^o x t])l^mj = /^on^mj- (44) 

Again we can neglect Eq. (44), because it coincides with (42). In the following we will not 
show third equations hke (41) and (44), because they are useless. Solution of all the equation 
is presented in Appendix B. Here for simplicity we limit ourselves only to a particular case of 
normal incidence of the waves. 

— The authors should not have allowed themselves 
to use the phrase "for simplicity". 
Referee of Am.J.Phys. 

4.2.3 A particular case of normal incidence 

In the case of normal incidence reflection and refraction are especially simple, because there 
is no splitting at reflection. We deflne geometry by three basis vectors n, I and t, where n 
denotes normal directed along z axis toward vacuum, while vectors I and t lie in the interface 
and deflne x and y axes respectively. The anisotropy vector a is supposed to lie in {x, z) plane 
at angle 6 with respect to n. 

A plane wave propagating along n {k = k/k = n) can have only two types of polarizations. 
It can be of a transverse mode with £i = ei = t and 'Hi = —nil, where ni = (see (2)), or 



it can be of mixed mode with £2 — ((16)), and 'K2 = n2{0)[n x a], where n2{0) = \Je2{9) 
(see (17)). 

Since there are no splitting at normal incidence the boundary conditions are simplifled. For 
the mode 1 (39) and (40) are reduced to 

[1 + 7^11] = ^ei, ni[l-yn] =^el, (45) 

and transmitted wave has fleld polarization E^i — t, £ei — I identical to that of the incident 
fleld. 

Solution of (45) is 

Pu = — T, ^el = — T, ni = y/€i. (46) 

rii + l ni + 1 

With these formulas we can immediately flnd reflection and transmission of a plane parallel plate 
of thickness D for an incident from vacuum electromagnetic wave with polarization — t: 

1 — Pi-^ exp{2ikiD) 1 — exp(2zfciL>) 
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where ki — koUi. We see that the incident wave with hnear polarization S^i — t parallel to 
that of mode 1 inside the plate does not change polarization after transmission through the 
plate. 

Now let's apply boundary conditions to the mode £2 incident normally on to the interface 
from inside the plate. Now the boundary conditions are reduced to 

{I • a)[l + ^22(^)1 = ^e2, k2{e){l ■a)[l- t22m = ^0 7^e2, (48) 

where factor (Z • a) appears because of projection of the polarization vector onto the interface. 
Therefore, since k2 — kon2{0), and n2{9) — \j^2{0) solution of (48) is 

^-(^^ = n2(^' ^^^^^^ = n2{e) + l ■ ^''^ 

From symmetry consideration we can immediately find refiection and transmission ampli- 
tudes for outside incident waves with unit polarization along I: 



"i^^M - ^ ^^-,^^^y (50) 

and therefore we can immediately find refiection and transmission of a plane parallel plate of 
thickness D for an incident from vacuum electromagnetic wave with polarization £2e — I- 

R2[0) - -T22W^_-^2^(^^^^p(2.^^(^^^^, 

T2{e) = eMiHO)D) . ^2\~/'i!!\ ^n,nV ^^^^ 

1 - P22{^) exp{2tk2{9)D) 

where k2{9) = kon2{9). We see that the incident wave with linear polarization £2e = I, which 
lies in the plane (n, a), does not change polarization after transmission through the plate. 

Now let's consider transmission through the plate of a plane wave ex.p{ikoZ — iuJot)£e with 
intermediate polarization: £e — at + pi, where lap + |^|^ = 1. The transmitted electrical part 
of the wave will be 

Etiz, t) = exp{iko{z - D) - iuot) [aTit + pT2{e)l] . (52) 

Since oTi = \aTi\ exp(i7i), and (3T2{9) = |/3T2(6')| exp(i72(^)), the real part of the wave (52) at 
some point z chosen for convenience so that ko{z — D) + 71 = 27rN with an integer N, looks 

Re{Et{z,t)) = \aTi\ cos(cuot)t + |/3T2(^)| sin(a;ot - ip)l, (53) 

where (p — ^2{9) — 71 — 7r/2. So, in this case the transmitted field has elliptical polarization. 

4.3 Reflection and refraction from outside an anisotropic medium 

Let's consider the case, when the half space at 2; < is vacuum, and that at 2; > is an 
anisotropic medium. The incident wave goes from the left in the vacuum. The wave function 
in the full space now looks 

*(r) = e(^ < 0)(e''^°'-?"^.+e^^°'- ^ "hp ^,^)+Q{z > 0) (e"^ '"-^ ^1^ + e'~^ '''^ 2'^^ 2j 

j'=e,m ' 

(54) 



\-n2{9) , 2 
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where j,j' denote e or m for TE- and TM-modes respectively, the term exp{ilc qv)'^ j with the 
wave vector ko — {k\\, ko± — ^k^ — /cjj) describes the plain wave incident on the interface from 
vacuum. In TE-modc factor = + contains — t and — [^o*]- In TM-mode 
factor "ipjn — '^m + '^m contains 

m 

— t and '^rn — — [^0^]- 
The reflected wave has the wave vector = {k\\, —kQ±), and fields ^g = t, ^g = [l^o*], 
— t, and = — [1« 0^]- The refracted field contains two wave modes with wave vectors 
fc 1 = {k\\,ki±), = ik\\,k2r±) and electric fields "^i = Ci = [aT^i] and '^2 = 62 = 
a — it 2{alt 2)^2{ d 2)1 ^A- Here n = k/k, ki± = y^ei/cg — kp and k2r± is given by (30). For 
incident TE-mode reflection pgg, pme and refraction r^g amplitudes (j = 1, 2) are found from 
boundary conditions 

{tt,)l^^, + {tt 2)l^2e = 1 + -fee, (55) 
(Z^l) T^le + (i^2) T^2e = -Ko±{l ' tee), (56) 
(t^l)^lg + (t^2) T^2e = fme, (57) 
(Z^l) T^lg + {it 2)^2e = -l^0±tme- (58) 

Solution of these equations is elementary and is given in Appendix C. 

4.4 Reflection and transmission amplitudes for a plane parallel plate 
of thickness L 

Now, whence understand what havens at interfaces, we can construct [17] expressions for 
reflection, R{L), and transmission, T{L), matrices for a whole anisotropic plane parallel plate 
of some thickness L, when the state of the incident wave is described by a general vector |"^o)- 
To do that let's denote the state of field of the modes Si and 62 incident from inside the plate 
onto the second interface at 2; = L by unknown 2-dimensional vector (109). If we were 
able to find \x) we could immediately write the state of transmitted field 

TiL)\-to)^r)^), (59) 
and the state of the field, refiected from the whole plate 

2(L)|to) = -^Ko) + r'E{L)n\lt), (60) 

where E{L), E{L) denote diagonal matrices 

^ feM^h.L) \ ^ feM^h.L) A 

y exp{ik2i±L) y ' \ exp(^fc2r±-t^) / 

which describe propagation of two modes between two interfaces. Here ki± — sjeikg — /cy , while 
k2r± and k2ii. are calculated according to (29) or (30) and (32), respectively. 

It is very easy to put down a self consistent equation for determination of |"^): 

1^) = E{L)T\^o) + E{L)nE{L)n\lt). (62) 

The first term at the right hand side describes the incident state transmitted through the first 
interface and propagated up to the second one. The second term describes contribution to the 



12 



Figure 5: Dependence of reflectivities |i?eeP and |i?m,eP of an anisotropic plate with ei = 1.6, 
T] = 0.8 and dimensionless thickness Lcu/c = 10 on angle (p of the plate rotation around its 
normal, when the anisotropy vector a is parallel to interfaces and at = is directed along 
The incidence angle 9 is given by sin ^ = 0.9. 

state of the itself. After reflection from the second interface this state propagates to 
the left up to the flrst interface, and after reflection from it propagates back to the point z = L. 
Two terms at the right hand side of (60) add together, which results to some new state. But 
we denoted it | af), and it explains derivation of the equation (60). 
Prom (60) we can directly find 

1^) =[i- E{L)nE{L)'R!]~'E{L)T\^o): (63) 
and substitution into (59) and (60) gives 

T{L) = i''' ]=r\i-E{L)nE{L)n] E{L)r, m 

\ me mm J 

= ^ + r'S(L)^[j - E{L)'^E{L)'S]~^E{L)T. (65) 
With these formulas we can easily calculate all the reflectivities and transmissivities 




IP 12 IP |2 \ / \T |2 IT 

I ^ee I I -'•'em | 1 j | -'■ ee | | -'■ em | 

I p 12 I p |2 I ' I It |2 IT" 

I J^Wne I ^mm / \ I me | | mm I 



(66) 



for arbitrary parameters, arbitrary incidence angles, arbitrary incident polarizations and ar- 
bitrary direction of the anisotropy vector a. In fig. 5 we present, for example, refiectivities of 
TE-mode wave from a plate of thickness L such, that Loj/c — 10. The anisotropy vector is 
parallel to interfaces. Therefore, its orientation with respect to wave vector fco of the incident 
wave varies with rotation of the plate by an angle around its normal. The transmissivities of 
the same plate in dependence on the angle (p are presented in fig. 6. 



5 D'yakonov surface waves (Dsw) 

Above we found that a surface wave can appear on the interface at total reflection of a mixed 
mode at some incident angles. This surface wave, however, is tied to the incident and reflected 
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Figure 6: Dependence of transmissivities |Teep and |TmeP of an anisotropic plate on angle of 
the plate rotation around its normal, with all the parameters the same as shown in caption of 
fig.5. 

wave and does not exist without them. D'yakonov in 1988 had discovered [18] (see also [18- 
21]) that on the surface of a uniaxial anisotropic medium there can exist free surface waves, 
analogous to elastic Rayleigh waves on a free surface. We will derive them with our tensor (2) 
and a little bit correct previous derivation by D'yakonov [18]. 

Let's again consider the space separated by a plane at z = to two halves, as shown in 
fig.2. The left part {z < 0) corresponds to anisotropic medium with dielectric permittivity (2) 
and the right part is an isotropic medium with dielectric permittivity e^. 

The surface wave is characterized by the wave function 

*(r, t) = [e(z < 0) (tAi exp(p,z) + -02 eMP2z)) + ©(^ > 0)iP, exp(-piz)] e"'\n-'''\ (67) 
where if} = S + 'H, 

Wi,2 = ^ ([Z X £1,2] - «gi,2[n X £1,2]) , ni = ^{[lyiSi]+iqi[nx£i]), (68) 

parameters pi,2,i provide exponential decay of the surface wave away from the interface. In (68) 
we also introduced dimensionless parameters qi^2,i — Pi,2,i/k\\. For £i^i we have 

qi,i = ^1 (69) 

respectively, where we denoted z = k^/k'^^. Parameters qi^i are positive reals when ei^iZ < 1. 
To find q2 for the field £2 we need to solve the equation 

+ r]{{l ■ a) - ix{n ■ a)f = zei{l + rj), (70) 

where x denotes q2, n is a unit vector of normal, directed toward isotropic medium, and I is 
a unit vector along fcy, which together with n constitutes the plane of incidence. From this 
equation it is seen, that §2 can be real only if vector a is perpendicular to n or to Z. In the 
first case the axis of anisotropy is parallel to the interface [18]: 

a — ail + Qft — cos 91 + sin 6t, (71) 
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and solution of (70) is 



?2 = ^1 + rjC0S^9 - €iz{l + rj). 
It is seen that 52(6') is positive real when 

l + V _ 



zei: 



= Z€2{9) < 1. 



In the second case 
and solution of (70) is 



1+7] cos^ 9 
a = ttnti + tttt = COS 4>n + sin 4>t, 



?2(0) 



COS^ 



(72) 

(73) 
(74) 
(75) 



Below we will show that in the second case free surface waves do not exist. 



5.1 Anisotropy axis is parallel to the interface 

When a is parallel to the interface, vectors £1,2, according to section (2), can be represented as 
C 

£1 = — — [a X k] = Ci (sin 6'n + iqi[a x n]) = Ci (sin^n + iq^ sin 91 — iqi cos 9t) , (76) 

Kh 



Co 



a — 



k{ak) 



l + V 



Co 



a - {I- iq2n){al)- ^ 

1-qi 



(77) 



k] l-ql + r]{l-af 
— C2 (iq2 cos 9n — q\ cos 91 -\- {1 — q\) sin 9t^ , 

where C2 = C2/{1 — qf), Ci,2 are some complex coefficients, and we used relation (72). 

In the basis n, I, t, shown in fig.2, polarization in isotropic medium can be represented 

as 

£i^ an + pi + -ft (78) 
with coordinates a, p and 7. Because of equation e^V ■ £i — 0, which is equivalent to 



iqia + 13 = 0, 



(79) 



vector (78) is reduced to 



£i = a{n - iqil) + -ft. (80) 
From continuity of t- and /-components of electric field at the interface we obtain two 



equations 



iCiqi sin 9 — C2q\ cos 9 = —iqia, 
—iCiqi cos 9 + C2{1 — q\) sin 9 = -f. 



(81) 



To get another two equations we need to use continuity of tangential components of magnetic 
fields. Substituting (76), (77) and (80) into (68) and neglecting common factor k\\/ko we obtain 

"Hi = Ci (^—iqi cos 9n + qf cos 91 — [1 — q^] sin 9t^ , 

I-L2 = ^"2(1 — ?i) (sin 9n + iq2 sin 91 — iq2 cos 9t) , 

Ui = (7n - iqal - a{l - q^)t) . (82) 
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Figure 7: Dependence v{6) = Cd{x)/c on angle between anisotropy axis a and direction k\\ 
of the surface wave propagation. The curve was calculated for ei = 1.6, r] = 0.4 and ej = 1. 



Continuity of I and t components gives 

qlCi cos 9 + iC2?2(l - ql) sin 9 = -iqi'y, (83) 

Ci sin 9 + iq2C2 cos 9 = cte, (84) 

where we denoted £ = ej/ei. If we exclude 7 and o; from these equations, we obtain a homoge- 
neous system of 2 equations for Ci,2 

qi cos 9{qi + gi)Ci + zC2((72 + gi)ei^sin^ = 0, 

iCi sin 9[qis + qi] - C2 cos 9 (sqj + qiq2j = 0. (85) 

The system of the two linear equations (85) has solution only if the determinant of its coefficients 
is equal to zero, which gives an equation for z = kQ/kfr. 



f{z) = qi{qi + qi) (eql + ^1^2) cos^ 9 - eiz{q2 + qi)[qi£ + qi] sin^ 6^ = 0. 



(86) 



Solution of this equation gives the speed of the Dsw cd{z) = C\fz. 

We derived equation (86) so scrupulously, to show that the result (86) slightly differs from 
the one presented by equation (8) of [18], and it is not reducible to Eq.(9) of [18], because 
solution of (86) exists even for < ei (ei is denoted e_L in [18]). More over it follows from 
(86) that the surface wave exists in much larger range of angles 9 than was obtained in [18]. 
For instance, in fig. 7 it is shown the dependence of ratio v{&) = Cd{z)/c on 0, it is seen that 
solution of (86) exists in the full range < ^ < 7r/2 for ei = 1.6, 77 = 0.4 and = 1. 

5.2 Anisotropy axis is perpendicular to the propagation direction 

When al-l then vectors £1^2, according to section (2), and "^1,2 according to (68) can be 
represented as 

C 

£1 = — — [a X fc] = Ci (sin 0n + iqi sin (f)l — cos 0t) , (87) 

k\\ 

82 = 02^ 1 — 5i + ?2 COS 0n + iq2 cos (pi + sin 9t^ , 
"Hi = Ci (— cos (pn — iqi cos 0Z — [1 — qf] sin (f)t^ , 
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"^2 — C2{1 — ql) (sin 0n + iq2 sin (f)l — cos 0t) , (89) 

and £i, Hi are the same as (80) and in (82) respectively. After performing the same procedure 
as above we obtain an equation for z — k^/k^i in the form 

fi{z) = (^^q2 + qi^ cos^ (p+{q2 + qi){i - qiqi) sin^ = 0, (90) 

which has no solution because all the terms in it are positive. Therefore the surface waves do 
not exist at such orientations of axis a, as was correctly pointed out in [18]. 

6 Conclusion-summary 

— / have no confidence that these authors can form 
professionally written, carefully researched work with 
a decent literature search. 
Referee of Am.J.Phys. 

In the case of uniaxial or biaxial anisotropic media we used for the tensor of dielectric 
permittivity Sij in the form 

£ij = ei [Sij + rjaiUj] , £ij = ei [Sij + rjaaiUj + rjbbibj] , (91) 

where ei is a parameter of isotropic part of the tensors, a, b are the unit vectors along axes 
of anisotropy, and r], r]a,h are respective anisotropy parameters. With such tensors we can 
immediately find for a plain wave £exp(ifcr' — iLot) with an arbitrary direction At = k/k of 
propagation analytical expressions for the polarization vector £ and wave number k{u:). In the 
case of uniaxial anisotropic medium we found that only two modes of linear polarizations can 
propagate inside it. One is transverse mode with 

£i = [a X /c], /ci = {u/c)^/e{, (92) 

and another one is the mixed mode (it has a component of polarization parallel to the wave 
vector) 

£2^a-K{K-a)'^, k2 = {u/c)^^), 62(g) = /^^^ ^ cosg=(a-/t). 
ei * l-\-ri COS"' i) 

(93) 

Next we considered reflection of obtained plain waves from an interface with an isotropic 
medium and had shown that reflection of every mode is accompanied by beam splitting, that 
the wave of mode 2 is in general reflected nonspecularly, and at some incident angles reflection 
of mode 2 can create a surface wave, which is bound to the incident and reflected waves of 
mode 2. The beam splitting at reflection can be spectacularly demonstrated with the help of 
light transmission through an anisotropic cone. 

After calculation of reflection from interfaces from inside and outside anisotropic medium we 
had shown an algorithm to calculate reflection and transmission of plane parallel plates without 
matching of the wave fleld at two interfaces. In the case of normal incidence on the plate of 
a plane wave with linear polarization transmitted wave in general has elliptical polarization. 
The form of ellipse changes with rotation of the birefringent plate around its normal and at 
two distinct orthogonal direction the ellipse reduces to linear polarization identical to that of 
the incident wave. 

Next we considered the free Dsw on the surface of anisotropic media. We corrected some 
error in derivation presented in [18] and have shown that Dsw exist in larger range of variation of 
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uniaxial crystal 



DSW 



Figure 8: Scheme of a possible experiment to observe Dsw. 



dielectric constants and in a larger range of angles between direction of surface wave propagation 
and direction of anisotropy vector a. 

To observe the D'yakonov surface waves it is possible to use the experimental scheme shown 
in fig. 8, which is different comparing to the one used in [22]. Disc of a uniaxial crystal with 
anisotropy axis a parallel to the surface can be pivoted around its axis to change angle between 
direction of Dsw propagation k = k\\/k\\ and the vector a. The Dsw is excited at frustrated 
total reflection in an anisotropic cone similar to that one shown in fig. 4 (here for simplicity 
we draw only one transmitted ray). Excitation takes place only when speed of the incident 
or reflected wave inside cone matches the speed of Dsw. Rotation of the cone around its axis 
permits some tuning of the speed. 

The second anisotropic cone identical to the first one detects Dsw, and the light transmitted 
into it through the small gap should be visible on a screen, as shown in fig. 8. 
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A Waves in a two axes anisotropic medium 

A two-axial anisotropic medium is characterized by two unit vectors a and b, and two anisotropy 
parameters e'^ and e'^. Therefore the tensor e has matrix elements 



In the last two equations we introduced notations rja = e'a/ei, and rjf, = e'f^/ei. For simplicity, we 
assume that a±b, introduce the orthonormal basis a, b, c = [a x b] and in this basis represent 



and equations (8), (9) and (10) take the respective forms 



(94) 



e£ = ei£ + e^a ■ £) + e;b(b ■ £), 

{k-£)+ r]a{Ki ■a){a-£)+ r]b{Ki ■ b)(b ■ S) = 0. 
[k'^ — k^tijE — A;^k(k ■ £) — kleir]aa{a ■ £) — kleir]hb(b ■ £) = 0. 



(95) 

(96) 
(97) 



£ = aa + /3b + ^yc 



(98) 
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with coordinates a, /?, 7, which are not completely independent, because (98) should satisfy 
(96). Substituting (98) into (96) gives 

a{k ■ a) + p{k ■ b) + ^{k ■ c)] + r]a{k ■ a)a + r]i,{k ■ b)p = 0. (99) 

Therefore 

7(/c • c) = -a{K • a)(l + 77J - /3{k ■b){l + rj^). (100) 
Prom (98) we also obtain that 

{£-a) = a, {£-b)=l3. (101) 

Now let's substitute (101) and (k • £) from (96) into (97) and multiply (97) consecutively by a 
and b. As a result we obtain a system of two linear equations 

{e[l + r]a{K-af]-eo^,{l + r]a))a + r]bk\K-a){K-b)f3 = (102) 

(k^[l + rib{K-by]-klei{l + rib))/3 + Vak\K-a){K-b)a = 0. 

The solution to this system exists if the determinant is equal to zero. This condition can be 
written as 

where 

ea,b(M= 1 "'^^^"^tl ' cos^„ = (Av-a), cose,^ {K-b). (104) 

1 + Va,b COS^ da^b 

The solution to (103) provides two different values of A;i^2, for which we find a, /3. After 
substituting the latter coordinates into (100), we obtain the last coordinate 7. Thus we find 
two different plain waves with wave vectors ki^2 — ^1,2'* and linear polarizations £1^2 (98). 



B Reflection from an interface from inside of the anisotropic 
medium 

Exclusion of T^ej from (39) and (40), and exclusion of T^mj from (42) and (43) gives two 
equations for "i^y, ~^2j-i which is convenient to represent in the matrix form 

Ul ■ %) + Ko±(t • "^1)^ ({I ■ %) + Ko±(t • "^2)^ \ ft, A = _ ■ ^^'^ + ''o^^* ■ '^^A 

(/^0±(t • %) - {I ■ tl)) (ACO±(t • %) - {I ■ ts)) I ~ W(* • %) - il ■ ^j)) ' 

(105) 

Solution of this equation is 

t2j) " ^ V - {'^0±{t ■ %) - {I ■ ti)) ((Z • %) + Ko±(t • ti)) ) U^(t • %) - {I ■ t,)} ' 

(106) 

Where D is determinant 

D^{{1- %) + Ko±{t ■ t,)) («o±(* • - {I ■ ^2)) - 

- ((Z • fc) + Ko±{t ■ ^2)) {Ko±{t ■ %) - {I ■ ti)) . (107) 
Substitution of these expressions into (39) and (42) gives refraction amplitudes e,mj 
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B.0.1 The most general case 

Above we considered the case when the incident wave has polarization vector Cj with unit 
amphtude. (We remind that vectors ej can be not normahzed to unity.) To find later refiec- 
tions from plane parallel plates we will need a more general case, when the incident wave has 
both modes with amplitudes Xi^2- To find amplitudes of refiected and transmitted waves in 
the general case it is convenient to represent the state of the incident wave in the form of 2 
dimensional vector 

then the states of refiected and refracted waves are also described by 2-dimensional vectors, 
which can be represented as 





$r) = ^M=7tV), l?o) = ^M=r'|i^), (110) 



where TZ and T" are two dimensional matrices 



We introduced the prime here and below to distinguish refraction and refiection from inside 
the medium and the similar matrices obtained for incident waves outside the medium. 

These formulas will be used later for calculation of refiection and transmission of plane 
parallel anisotropic plates. In the case of a plate we have two interfaces, therefore we need also 
refiection and refraction at the left interface from inside the plate. They can be easily found 
from symmetry considerations. Their representation is obtained from (106) — (108) by reverse 
of arrows and change of the sign before ko±. After this action we find 

^' = ( h t'^], r'^i Z'' ] ■ (112) 

V T21 ¥22 J \ rmi rm2 J 

Refiection from outside the medium is to be considered separately. 

B.0.2 Energy conservation 

It is always necessary to control correctness of the obtained formulas. One of the best controls 
is the test of energy conservation. One should always check whether the energy density fiux 
of incident wave along the normal to interface is equal to the sum of energy density fiuxes of 
refiected and refracted waves, and the most important in such tests is the correct definition of 
the energy fiuxes. In isotropic media it is possible to define energy fiux along a vector n as 

(fc-n) c e£^ + 
(n.[5x^])_ 

In isotropic media both definitions are equivalent, because l-L = [k x £]/ko, and {k ■ £) = 0. 
The first definition looks even more preferable since the second one can be written even for 
stationary fields, where there are no energy fiux. 

In anisotropic media only the second definition is valid, and because in mode 2 the field £ 
is not orthogonal to k, the direction of the energy density fiux is determined not only by wave 
vector, but also by direction of the field £ itself. 
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C Reflection from an interface from outside of the anisotropic 
medium 



Exclusion of "j^gg and ~^ j^e leads to 



%±{Ml) + {it I 



(Ml) 



Ko±{M2) + {it 



{M2) 



and the solution 



^ _ {ko±{M2) + {it 2)) - {Ko±{tt2) - {1% 
;t2e) ~ D,\- (ko±{M,) + {it,)) («o±(*^l) - (i^l)) 



where Dp is determinant 



De - {l^O±{tti) - {Ml)) {no±{M2) + {it 
- {l^0±{tt2) - {M2)) {ko^{Mi) + {it I 



Substitution of l^je into (55) and (57) gives 



.1 



ee 



{tti) {tt2) 

{Ml) {M2) 



.0; 



In the case of incident TM-mode we have boundary conditions 

{Mi)l^lm + {M2)^2m = 1 + -f mm. 

{lti)'tim + {lt2)l^2m = «0±(1 " "^mm)> 
{tti)l^lm + (t"^2)^2m = 'fern. 
{l'Hi)l^lm + {l'H.2)^2m = l^0±l^ em.- 

Exclusion of T^^e ^^'^ 1^ mm Icads to 



>0±(*^l) 



{itl) 



[ko^{M2) 



{it 



Ko±{tti)-{Mi)) (Ko±{tt2)-{M 



Therefore 



Im 



Ko±{M2) 



~T^lm 
J^2mi 



{it. 



where Dr, 



^«:o±(t^2) - {1% 

Dm \ -\Ko±{tti) - {Ml)) (ko±{Mi) + {it I 

-De (117). Substitution of l^jm into (119) and (121) gives 

/ {tti) itto) \ ^^^m\ fO 



{Ml] 



{tc.2) 

{M2) 



2koj 
. 



. , 



, , 



2ko± 
, , 



(115) 



(116) 



(117) 
(118) 

(119) 

(120) 
(121) 

(122) 

(123) 



(124) 



(125) 



In the general case, when the incident wave has an amplitude ^ f. in TE-mode and amplitude 
in TM-mode, the state of the incident wave can be described by two-dimensional vector 



I to) 



m) 



(126) 
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and the states of reflected and transmitted waves can be represented as 



f „> = (I J = ^l7«>. it) = (I j = rit„>, (127) 



where 7?- and 7" are the two dimensional matrices 



1^' V r = f J'" J" V (128) 

V Pme Pmm J \ T 2e r 2m J 
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